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Abstract 

We discuss the effects of resonant z/g —>■ Ug and Uf. —» z/^ (z/^ is a sterile 
neutrino) conversions in the dense medium of a supernova. In particular, we 
assume the sterile neutrino Ug to be in the hot dark matter few eV mass range. 
The implications of such a scenario for the supernova shock re-heating, the 
detected Ue signal from SN1987A and for the r-process nucleosynthesis hypoth- 
esis are analysed in some detail. The resulting constraints on mixing and mass 
difference for the z/g — Us system are derived. There is also an allowed region 
in the neutrino parameter space for which the r-process nucleosynthesis can be 
enhanced. 
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1 Introduction 

The possibility that hght sterile neutrinos can be mixed in the leptonic charged current 
seems to be the very appealing from the point of view of the present anomalies observed 
in the neutrino sector: the solar [|l| and atmospheric [Q] neutrino problems as well as the 
need for a few eV mass neutrino as the hot dark matter in the Universe 0, ^ . Barring the 
possibility that the three active neutrinos are nearly degenerate in mass ||^, the simplest 
way to simultaneously account for these observations is to postulate [] the existence of a 
light sterile neutrino 0, ^ 0. Moreover, some of these scenarios may also account for 
neutrino oscillations between the electron neutrino and muon neutrino, as possibly hinted 



at the LSND experiment ||T0| 



The conversions between neutrino species may change significantly the phenomena 
occurring in supernova. Of particular interest is the conversion to a sterile neutrino {SU{2)l 
singlet), since the sterile state does not interact at all with the supernova matter. 

Here we focus on the resonant conversions of electron neutrinos (or anti-neutrinos) 
to sterile neutrinos outside the neutrinosphere. For mass difference Sm? < lO^eV^ the 
MSW resonance will occur in these regions of the supernova where neutrinos freely stream. 
Note that for this mass range the conversions to sterile neutrinos in the inner core can 



be neglected, for all values of the mixing angles |TT[]. Hence, in the absence of other 
non-standard interactions, such as a large magnetic moment, the sterile neutrinos are not 
emitted from the supernova core. We will work always under this assumption. Moreover, 
we assume that the mass eigenstate consisting mainly of the singlet state is heavier, with 
mass between 1 eV and 100 eV. The opposite case would require however special forms of 
the mass matrices in order to avoid the constraints from neutrino-less double beta decay. 
This possibility certainly exists but, for definiteness we do not consider. Finally, we assume 
that no other neutrino conversions, such as flavour conversions, take place. Generalization 
to the case of three- neutrino flavours is straightforward but model-dependent. 

The effects of such active to sterile transitions have been previously discussed in 



1^ , p!3| , |IJ], mainly concentrating on the effects for neutrino observations at terrestrial 
detectors. In the present paper we consider the implications of resonant u^. -^ Ug or z/g -^ i?s 
conversions in the dense supernova medium on the neutrino re-heating of the shock wave 
on the Ug signal in terrestrial detectors, as well as on supernova heavy-element 



nucleosynthesis |]T7 



The effect of the electron neutrino to sterile neutrino conversion on supernova nu- 



^The measurement of the Z boson width at LEP has provided a hmit on the number of active neutrino 
types with mass smaller than Mz/^, Ni, — 3.09 ±0.13 H. However, it cannot preclude the existence of 
light sterile neutrinos. 



cleosynthesis (via rapid neutron capture process or so-called r-process) has not yet been 
discussed with sufficient detail. So far most studies have concentrated on the case of active 



neutrinos, giving stringent limits on electron, muon and tau neutrino conversions |ll8| , [19 
Nevertheless, it has been suggested [^ that a resonant conversion from muon neutrinos 
to sterile neutrinos, with a subsequent conversion between electron and muon neutrinos, 
would enhance r-process nucleosynthesis. The required conversion pattern appears natu- 
rally in some specific models 0. 

Neutrino conversions may also influence the explosion mechanisms of the supernova 
in different ways. Although the main interest lays in the transitions among the active 
neutrino species [B^, it has also been suggested that conversions between sterile and active 



neutrinos may enhance the explosion |2^ . However, for the mass range we are considering, 
such an enhancement is not possible, the conversions weaken the shock wave, preventing 
the explosion. 

In Sec. 2 we give a quick reminder on the picture of the neutrino propagation in matter 



and on the resonant v^ — v^ conversion [23|. In Sec. 3 we give a qualitative discussion of 



the electron concentration Y^. in supernova and on the effects of the active-sterile neutrino 
conversions on Ye and, in turn, on the neutrino evolution itself. This can lead to non-trivial 



feedback effects [0, ^, ^. Sec. 4.1 discusses the implications of z/g — z/^ conversions for the 
neutrino re-heating mechanism. In Sec. 4.2 we analyse the impact of our scenario in the 
later epoch of supernova evolution (few seconds after the core bounce) for the supernova 
(anti)neutrino detection (Sec. 4.2.1) as well as for r-process nucleosynthesis (Sec. 4.2.2). In 
Sec. 5 we summarize our results and discuss their significance, by comparing the supernova 
limits we derive with the laboratory and nucleosynthesis limits on z/g — Vg conversions. 



2 The active- sterile neutrino resonant conversion 

In our discussion we only consider the conversion channels z/g — > z/^ and z/g — * z/^ where 
Vs {vs) is the sterile neutrino ^. For the sake of simplicity we will not consider the effect of 
^e ^^ v^i,T and z/g ^^ z/^^T- conversions. In the following we consider the 5m? = m\ — m\ > 
case, corresponding (for sufficiently small mixing angle) to the situation in which the heavier 
state is mostly the sterile neutrino. 

The evolution of the z/g — Us system in the matter background is determined by the 
Schroedinger equation 

A(^.] = ("' »-](■'']. (1) 



dr \ z/, / \ iJg, H. 



es 



z/.. 



^In the ultra-relativistic limit Vg and i^g have opposite hehcity. 



He = Ve-—^cos29, H,= V, + --—cos29, 

Stti 
H,s = ^;^sin2^, (2) 



where the effective potential Ve for Ue arises from the coherent forward neutrino scattering 
off-matter constituents [^ and is given by Q 

K . ^(n-iv;,^^,3>;-i), (3, 

J- e ^ I ; -' n ^ ^ e ■ 

Tie + Tin 

Here Gp is the Fermi constant, p is the matter density, rriM is the nucleon mass and Ue 
and Tin are the net electron and the neutron number densities in matter, respectively. Note 
that charge neutrality Yp = Yeis assumed and that there is no potential for z/^, i.e., ^ = 0. 
For the system z/g — > Vg the matter potentials just change their sign. 

The resonance condition reads as: 

K = ^cos2e. (4) 

Let us remind that for 5m? > 0, either the conversion z/g — > Ug (for Ve > i.e. 
Ye > 1/3) or z/e — > z/g (for VJ= < i.e. Yg < 1/3) takes place. This is important because, as 
we will see later, in the region above the neutrinosphere the matter potential Ve changes its 
sign due to the different chemical content. For our later discussion, it is instructive to know 
the profiles of the matter density and of the electron fraction Ye outside the neutrinosphere. 
In Fig. la and lb we plot these quantities for t < 1 s post-bounce (in short p.b.) and 
t > 1 s p.b. as given by the Wilson supernova model |^. We can see that the electron 



concentration Yg is rather low just near the neutrinosphere. Ye ~ 0.1 and 0.01 for the 
earlier and later epoch, respectively and far away it increases to values > 0.4. On the 
other hand, the matter density p exhibits a monotonically decreasing behaviour. In Fig. 
2a and 2b we plot the modulus of the effective matter potential Ve using the matter density 
and Ye profiles as given in Fig. la and lb. The position where Ye takes the value 1/3 (i.e. 
V^ = 0) is indicated by r*. This position corresponds to r* ^ 160 km and 12 km for the 
earlier and later epochs, respectively. Clearly, the effective potential Ve changes its sign 
from negative to positive at the point r*. 

The resonance condition in the eq. (^) provides the Sm^ value for which neutrinos 

with some given energy can experience the resonance for a certain value of the potential 

■^The effective potential should also contain contributions from the neutrino background. We have 
ignored them since the neutrino densities in the relevant regions are at least one order of magnitude 
smaller than the corresponding electron densities, and the neutrino terms in the effective potential involve 
an additional suppression factor because most neutrinos travel almost in the same direction. 



(or equivalently, at some position r). For the sake of convenience, in the right ordinate 
of Fig. 2 we have also indicated such corresponding values of Sm"^ for typical neutrino 
energy E = 10 MeV. We see that for 6m'^ > 10^ eV^ only Ue — > i^s conversions can take 
place, and this happens in the region where Yg < 1/3. On the other hand for smaller 
values, dm"^ < 10^ eV^ there can occur three resonances. The z/g's are first converted, say 
at ri < r*, then there are two resonance points at r2 and r^, where r^ > r2 > r* (i.e. in the 
region where Y^ > 1/3), for the z/g <-> Ug channel. In order to illustrate this more explicitly 
we plot in Fig. 3 the schematic level crossing diagram for z/g — Ug and Ug — Us system, 
assuming the mixing angle to be small. From this figure, it is also manifest that the z/^'s 
originated from the first z/g conversion (at r2) can be re-converted into z^g's at the second 
resonance (at r^). In our subsequent discussion, we will employ the simple Landau-Zener 



approximation |^, ^ to estimate the survival probability after the neutrinos cross the 
resonance. Under this approximation, the z/g (or z/g) survival probability is given by (in the 
case of small mixing angle) 



exp 



exp 



2 Li 



-2 X 10"^ X sin^2^ 






/lOMeVy 
V E, ) 



dVp km 



\/ 

dr eV 



(5) 



where L^*^ is the neutrino oscillation length at resonance. Notice that for 5r /U^ > 1 the 
resonant neutrino conversion will be adiabatic [^. We can expect the maximal sensitivity 
to the mixing angle for Sm^ = lO^eV^. From Fig. 2a and 2b one can estimate that 
in this case the resonance occurs at high density p ~ 10^^ — lO^^g cm~^ just above the 



neutrinosphere, where the gradient 



dVe 

dr 



1 — 10 ^eV km ^. From the eq. 



we can 



estimate that the conversion will be adiabatic for sin^ 20 > 10 ^ — 10 ^. 

Ug channel, the z/g survival probability 



Due to the double resonances for the z/g 
after the second resonance is given by 



PU 



P{r2)P{rs) + [l-P{r,)][l-P{r,)], 



(6) 



where P{r2) and P{r-i) are the survival probabilities calculated according to the eq. (^ at 
r2 and rs, respectively. 



3 The feedback induced by the neutrino conversion 



All the neutrino species emitted from the neutrinosphere have approximately the same 
luminosity L^ after a few ms p.b., characterised by a thermal Fermi distribution with 
temperature T^, and zero chemical potential. The typical duration of the neutrino emission 
is about 10 seconds. 



3.1 Neutrino emission and absorption reactions and Y^ profile in 
supernova 

In the region outside the neutrinosphere, due to the intense neutrino radiation, the electron 
fraction is determined by the neutrino capture by nucleons and by their reverse processes: 

Ue + n ^ p + e~ , (7) 

Ue + P ^^ n + e~^ . (8) 

In particular Yg above the neutrinosphere is approximately given by [l 

Y ~ Ae+„ + X„^n /„x 

X-I-X-I-X--I-X' ^^ 

The neutrino capture rates depend essentially on the neutrino luminosity Ly and energy 

Kn ~ / (t)\E,)a,N{E,)dE, ex -j^{El) oc L,{E,) , (10) 

where (z/, N) = (z/g, n) or (z/g, p)- The neutrino luminosity is given by the black-body 
surface emission formula, L^ = 5.6 x lO'^^r^ T^ ergs/s (here r,y is given in km and T,^ in 
MeV). For our purpose it is sufficiently accurate to assume the neutrino differential flux 
(jP{Eij) to be given by 

^ ^ "^ 47rr2/E3dE,/(e^-^/^-' + l) " ^ ' 

On the other hand, the rates for the inverse processes depend strongly on the matter 
temperature T and are given as 

TT^ Jo exp[(Ee =F MlT\ + 1 

where a^N is the electron or positron capture cross section with (e, A^) = (e~, p) or (e^, n). 
(In the above formula the negative sign in front of /ig is for electrons and the positive sign 
is for positrons.) 

In Fig 4a and 4b we have plotted the above rates X^n-, Xsn (for the earlier and later 
epochs, respectively) as a function of the distance from the centre of the star. For conve- 
nience we have also plotted the matter temperature profile by the solid line (see the right 
ordinate scale). One can see from Fig. 4a and 4b that close to the neutrinosphere the rate 
Ae-p dominates over the others, and hence Y^ is less than 1/3 (see the eq. (^), as can 
be seen from Fig. la and lb. As r increases Ag-p decreases faster than Aj^tv and hence Y^. 
increases. At some point r r^ r* where Ag-p ~ X,yN, the fraction Ye takes the value 1/3. In 
the region r > r* where the neutrino absorption reaction is high enough to dominate the 



proton-to-neutron ratio, the electron fraction is given by 



y; ~ ^"-"^ ~ \ . (13) 



assuming the fluxes to be given as in the eq. (]Tl[). For the typical energies {Ey^) ~ 11 MeV 



and (-Epe) ~ 16 MeV at t > 1 s p.b., we find Ye = 0.41, in good agreement with numerical 



supernova models |T^. The above discussion explains well the behaviour of the Y^ profiles 



shown in Fig. la and lb. 

3.2 Feedback effect on neutrino conversion 

Now let us come to the main point concerning the effect of the neutrino conversion upon 
the electron content in the matter. Suppose that at some position tq and time to neutrino 
conversions z/g -^ Vs or Ve —* ^s occur. From the eq. @) and eq. (|^) we see that the 
conversion probability depends on the value of Yg and its derivative at the resonance 
position tq. The neutrino conversions Vf. -^ Vg or Pg -^ Vg reduce the value of Xy^n or Ap^p 
and this could lead to the modification of Y^ in the region r > tq (see the eq. (^). If 
the modification of Ye is fast enough it could affect any subsequent neutrino conversion 
occurring at some position r > tq and time t > Iq. This subsequent neutrino conversion 
could again modify the value of Ye-, and so on. In this way, neutrino conversion rates and 
the electron fraction continuously affect each other during the neutrino emission, leading 
to non-trivial feedback phenomena for the neutrino transitions. These should in principle 
be taken into account self-consistently in neutrino conversion studies. It is well known that 
in the early universe feedback effects on the neutrino background should be reflected in the 



evolution of neutrinos [24] and such effects have been studied in ref. [25|. On the other 
hand, in the context of supernova physics, feedback effects on the neutrino background 
have recently been studied in ref. |2^. In what follows we will consider the analogous 
effect for the electron abundance Yg- 

From our discussion in Sec. 3.1, we can immediately understand that the feedback 
is operative only when X^^ ^ X^n (see the eq. (||)). Right after the bounce (t < 1 s p.b.) 
the relevant resonance layers lie so far away from the stellar core that the Aj^tv is too small 
since 0° oc rl/r"^ (cfr. eq. (0) ) and as a result feedback effects would be small. On the 
other hand, later on (t > 1 s p.b.) the neutrinosphere shrinks so much that Aj^at is high at 
the relevant resonance positions (now they lie much closer to the centre) and, as a result, 
the feedback effects would be potentially important. 

The neutrino conversion could be affected by the feedback either by shifting of the 
resonance position, or through a change in the adiabaticity. Let us first see the effect on 
the resonance position. For simplicity let us assume that Y'g is determined only by X^n as 
in eq. (|T3|). In the region where Y'g > 1/3 (Y'g < 1/3) the z/g — > Vg {pe -^ i^s) conversions 

6 



occur and this tends to decrease (increase) Y^ in the region above the resonance layer 
because X^^n (-^Pep) decreases. This could affect the resonance condition for the subsequent 
neutrinos, leading to a shift of the resonance position. Moreover, if the value of Y^ becomes 
less or equal to 1/3, the subsequent z/g -^ ^s conversion would be suppressed because the 
resonance condition is not satisfied. This is therefore a negative feedback. Analogously, a 
negative feedback is also present for the z/g —^ ^s channel. 

Let us now turn to the feedback effect on the adiabaticity. The neutrino conversion 
would change the potential gradient gradually, as neutrinos of different parts of the energy 
spectra reach the resonance. In the case of Ve -^ ^s conversion, the depletion of the anti- 
neutrino fiux drives Y^ up by decreasing the anti-neutrino capture rate, and so it tends 
to raise the absolute value of the gradient of the potential. This weakens the adiabaticity 
of the neutrino conversion. As for the z/g -^ Vg conversion, in the increasing part of the 
potential in Fig. 2, due to the decrease of Y^ caused by Ve to Vs conversion, the shape of 
the potential would be fiattened, leading to better adiabaticity. On the other hand, for 
the decreasing part of the potential (see again Fig. 2), the potential gradient would be 
steepened leading to worse adiabaticity. 

Let us now consider the relevance of the feedback in terms of the 5m^ involved. 
We can distinguish three different ranges of drn^: (i) 10^ eV^ < dirP' < 10^ eV^, (ii) 10^ 
eV^ < 5m^ < 10^ eV^ and (iii) drn^ < 10^ eV^. For the mass range (i) and (ii) (see Fig. 
2), only anti-neutrino conversions take place, while for the range (iii) both neutrino and 
anti- neutrino conversions take place. In range (i), the transitions occur close enough to the 
neutrino sphere, where the e~p -^ z/gU reaction dominates over the corresponding neutrino 
absorption reactions. In this case Y^ around the resonance position is not affected by the 
'^e —^ ^s conversion. Hence the feedback is irrelevant for this mass range. 

On the other hand, for the range (ii), the conversions occur farther away from the 
neutrinosphere where the e~p -^ Ueii rate becomes comparable to that of u^p — > e+n, 
especially in the later epoch. Hence the anti-neutrino conversion would be somewhat more 
affected by the feedback than in the previous case. However, the effect is still not large, 
since Ag-p is not small compared to Aj^at. 

Finally, for the range (iii) the situation is more complicated. In this case three 
resonant conversions can occur: one in the z/g — > z/g channel and two in the z/g ^-» z/g channel 
(see Fig. 3). For the z/g -^ Vs resonant conversion (at ri < r*) the preceding conclusion 
still holds and the feedback is not important. Similarly the next conversion z/g -^ Ug (at 
r2 > r*) is not too affected by the feedback due to its proximity to the core, as before. 
Moreover, this z/g -^ z/^ conversion is expected to have the same degree of adiabaticity as 
that of the previous anti-neutrino conversion, due to roughly the same steepness of the 
matter potential. As a result P{ue ~^ ^e\ ^i) ~ P{ve -^ ^e! ^2) and, to some extent the 



change in Y^. at r2 will be compensated by that at ri, leading to a small net feedback 0. 

Now let us come to the last resonance in the channel Ue ^^ i^s, at r^. This conversion 
could be significantly affected by the feedback because in the region around r^ (:» r*) Y^. is 
determined only by the X^^n reaction rates, since the inverse rates are small, as seen from 
Fig. 4b. For the sake of discussion, let us first consider the case where the conversions 
at ri and r2 are sufficiently adiabatic. In this case the resonance at r^ is also expected 
to be very adiabatic, since the slope of the potential around r^ is much less steep than 
that for the previous resonances. Therefore, the u^. fiux would be completely recovered by 
the last conversion Ug —>■ Ue at r^. This last conversion Ug —>■ Ue would remain adiabatic 
because the increase of Yg due to the conversion would not violate the adiabaticity but 
simply shift the resonance position. As a result we conclude that it would not be affected 
by the feedback. In the opposite case, when at ri and r2 the conversions are not adiabatic 
(P ~ 1), the feedback effect, due to the last u^ <-> z/^ transitions, could be very important. 
The depletion of the electron neutrino flux and the unchanged Ue flux would lower the value 
of Ye and hence Ve, with the effect of shifting the resonance position r^ at lower values. At 
the same time the potential gradient becomes steeper and hence the conversions become 
less adiabatic, leading to the suppression of the z/g — >■ Ug transition. In this case we would 
have a net feedback effect on the neutrino conversion. This situation occurs for rather 
small mixing angle sin^ 26' < 10~^ and for 6m'^ < 10^ eV^. We can have a quantitative 
insight of the signiflcance of the effect from the eq. (|13D by simply requiring Y^ > 1/3. For 
simplicity, let us assume Pp^(ri) ~ PuX''^^) ~ 1- Thus we can deduce a lower bound for the 
survival probability at rs, P^^(r3) > 0.7. In other words, the second z/g — > Vs conversion 
would be stopped when about 30 % of the lower part of z/g spectra are converted to z/^ 
(because the lower energy neutrinos undergo resonance before the higher ones) even if the 
conversion is initially very adiabatic. 

Instead of taking into account this feedback effect upon the neutrino evolution, we 
have simply stop the z/g — ^ z/^ conversion when Yf, reaches the value of 1/3. 

We conclude that the feedback effects would be small and could be neglected in the 
earlier epoch in the region relevant for our discussion. As for the later epoch, the feedback 
may be more relevant especially for the range dm? < 10^ eV^. 



''This observation is strictly true in the case of absence of the feedback. However, since the feedback 
effect may sUghtly decrease the adiabaticity of anti-neutrinos and increase the adiabaticity of neutrinos, 
the conversion probabihty for neutrinos may be a httlc larger than for anti-neutrinos. 



4 Constraining Neutrino Parameters 

In this section we are mainly concerned with the iniphcations of active-sterile neutrino 
conversions for supernova physics. For the earlier epoch of supernova evolution active- 
sterile neutrino conversions would suppress the shock re-heating. For the later epoch 
active-sterile neutrino conversions could suppress the detected z/g signal from SN1987A. 
On this basis we derive stringent constraints on the neutrino parameters. On the other 
hand, these conversions could affect the r-process nucleosynthesis in the later epoch, either 
to suppress it or to enhance. In the first case we again analyse the restrictions on neutrino 
parameters. 

4.1 Earlier Epoch: t < 1 second after core bounce 

In the following we consider only the epoch after the core bounce and the neutrino evolution 
in the regions outside its neutrinosphere. Indeed, for the range 6m'^ < 10^ eV^ neutrino 
transitions in the dense matter of the core {p ~ 10^^ g/cm^) are strongly suppressed |ll|, |13 



At this time a reflected shock wave is formed between the inner homologously col- 
lapsing part and the supersonically falling outer portion of the initial iron core. This shock 
propagates through the mantle and may result in mass ejection when it reaches the sur- 
face of the star. In fact the shock suffers energy loss due in particular to the emission of 
neutrinos which prevents a successful explosion of the star. As soon as the shock wave has 
passed through the neutrinosphere there is a large burst of UeS. Subsequently on a time 
scale of several seconds after the core bounce the emission of all neutrino species drives 
the evolution of the star to the final cool and neutron star. The epoch within 1 second 
after the core bounce is rather important for the re-heating of the shock wave. In the 



delayed explosion mechanism |15, 16] the neutrino energy deposition, occurring between 
the neutrino sphere and the site where the shock is stalled, can re-start the shock and 
power the explosion. As can be seen from Fig. 4a, at r ~ 300 -v- 400 km away from the 
neutrinosphere the neutrino absorption rate on free nucleons dominates over the capture 
of electrons. At this position the energy transfer from neutrinos to the matter takes place 
so as to help the shock Q 

In the absence of neutrino conversions the corresponding energy gain (per nucleon) 

f ^ m is 

m « ^ (14) 

where L^, is the total u^ + i^e luminosity and cr^^ ~ 9 x 10^'*'^ x E^ cm^. For Lj, ~ 3 x 10^^ 



^ Other mechanisms for energy deposition, hkc neutrino scattering off-electrons or neutrino-anti- neutrino 
annihilation are less efficient |pQ] and thereby we neglect them in the following. 



ergs/s, E^ = 10 MeV and r = 300 km, we find E{t) ~ 20 MeV s~^. This rate seems 
to be large enough on the time scale of 0.1-0.2 second if compared with the gravitational 
potential (per unit mass) GNMr/r ~ 7-10 MeV of the material stopped behind the shock 
(Mr ~ 1.5Mq is the included mass). Thus the neutrino energy transfer can help the 
material to overcome the gravitation of the star and so to escape it. The success of this 
scenario strongly depends on the neutrino luminosity. Clearly the sterile conversion z/g -^ Vs 
or z/g — * Vs would imply a depletion of active neutrino luminosity and thereby can spoil the 
re-heating process |1^ . 



We have calculated the ratio R of the neutrino heating rate in the presence of Ve -^ ^s 
and Ve —^ ^s transitions to the corresponding rate in the absence of such transitions. 
Following ref. [^ we use an approximate expression for i?, in which we neglect the re- 



emission of neutrinos by the heated matter, leading to 

EMt) ~ j E,a,M(t)\E,)dE, , El^{t) r^ j E,a,MP{E,)(t)\E,)dE, . (16) 

where the primed quantities Y' and Y^ stand for the proton and neutron abundances 
calculated in the presence of active-sterile neutrino conversions. 

In Fig. 5 we plot the iso-contour for different values of the ratio R in the parameter 
space {5m?, sin^ 29). Requiring a moderate effect R > 0.9, one can exclude sin^ 26 > 10~^ 
for Sm"^ ~ 10^ eV^, whereas for smaller mass, 6ir? 1 ^ 10 eV^, the bound on the mixing is 
weaker and lies in the range sin^ 29 > 7 x 10^^ ^ 5 x 10^^. Note that our bounds are in 



qualitative agreement, though slightly more stringent than those found e.g. in ref. |13 . 



4.2 Later Epoch: t > 1 second after core bounce 

For the later epoch we consider the effect of active-sterile neutrino conversions both on the 
z/e signal as well as the r-process nucleosynthesis and analyse the possible restrictions on 
neutrino parameters. 

4.2.1 Implications for the detection of SN1987A Ug signal 

The Kamiokande II and 1MB detectors observed 11 and 8 Ug events, respectively, from 
SN1987A pT|, R21. This is in agreement with the theoretical expectations, which predict 



that almost all of the released gravitational energy is radiated in all neutrino and anti- 
neutrino flavours. Signiflcant conversion of z/g's into a sterile neutrino would be in conflict 
with this evidence. We can just apply this consideration to constrain the neutrino mixing 
and mass difference. 

10 



We plot in Fig. 6 three contours of the z/g survival probability P for the Ve —>■ i^s 
conversion, in the (5m^,sin^2^) parameter space. The upper line is for P = 0.1, the lower 
one is for P = 0.7 and that in the middle corresponds to P = 0.5. If we assume that the 
successful observation of the SN1987A signal implies that at least 50 % of the expected 
z/e signal has been detected, one can conclude that all the portion above the contour of 
P = 0.5 is ruled out. For 6171"^ = 10^ eV^, the range sin^ 26 > 5 x 10~^ is excluded whereas 
for Srn^ = 1^10 eV^ the non-adiabatic character of the conversion implies a much looser 
bound, sin^ 29 > 10~^. Our results are again in qualitative agreement with those in ||12l 113 



Now we would like just to briefly comment on the relevance of the large volume 
detectors, such as Superkamiokande |^ or SNO |^| , aimed to detect supernova neutrinos. 



A galactic supernova event would produce e.g. about 5000 events through u^p reaction in 
Superkamiokande detector. Such a huge statistics may not only allow to determine the 
neutrino flux with good accuracy but also may provide the necessary sensitivity to measure 
e.g. the neutrino energy spectrum. The resonant conversion between electron neutrinos 
and sterile neutrinos may show up as a deficit of neutrino events, distortion of spectra, 
time dependent flux, etc. The absence of a deficit in the expected number of events can be 
used to further constrain the neutrino parameters. In Superkamiokande detector a number 
of z/e events not larger than 2500 would disfavour all the region above the iso-contour of 
P = 0.5 in Fig. 6. However the observation of a possible deficit on electron neutrinos or 
anti-neutrinos is not in itself a distinguishable signal for any specific conversion mechanism. 
Moreover, the precision of any conclusions that can be drawn from these considerations is 
limited by the uncertainties in the theoretical neutrino fiirxes. 

A better signal is the distortion of the neutrino energy spectrum that would arise 
from the neutrino conversion. The energy dependence of the adiabaticity condition may 
already cause a mild distortion, typically neutrinos of lower energies have larger conversion 
probabilities. That is, however, a general feature of the resonant conversion irrespective 
of the specific neutrino channel. However, we can envisage a possible distortion caused by 
the feedback effect. As discussed in Sec. 3, the feedback effect may cause a 30 % reduction 
of neutrino events in the low-energy portion of the neutrino spectrum (with a sharp cutoff 
energy close to 10 MeV). In other words, there will be a tendency for neutrinos above this 
energy to be blocked from converting. In order for this to happen, the following conditions 
are to be satisfied: 1) the first resonances (for electron neutrinos and anti-neutrinos) must 
be non-adiabatic i.e. for Sm"^ < 10^ eV^, sin^ 29 < 10~^; 2) the last resonance should be in 
the region where the electron fraction is determined by the neutrino absorption reactions, 
i.e. for Sm? > 1 eV^; 3) the last resonance should be adiabatic. A more detailed and 
careful analysis will be given elsewhere [p5| . 
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4.2.2 Implications for r-process Nucleosynthesis 

The implications of resonant neutrino conversions into active neutrinos for the supernova 



nucleosynthesis have been recently investigated in a number of papers [0, M 



The r-process is responsible for synthesising about half of the heavy elements with 
mass number A > 70 in nature. It has been proposed that the r-process occurs in the region 
above the neutrinosphere in supernovae when significant neutrino fluxes are still coming 
from the neutron star |T^. A necessary condition required for the r-process is Y^ < 0.5 



in the nucleosynthesis region. The Yg value at large radii above the neutrinosphere, where 
the r-process nucleosynthesis takes place, is determined, as we have discussed in Sec. 3, 
only by the neutrino absorption rates Xi^^n and Xp^p- Therefore, Y^ in the nucleosynthesis 
region is approximately given by the eq. (0). We have also learnt that the presence of 
neutrino conversions into a sterile state can affect the corresponding Yg. Thereby, in the 
nucleosynthesis region we can write Yg as follows 

As we already noted the u^. —>■ Vg conversion leads to a reduction of v^ luminosity and hence 
to an increase of Y^, whereas the z/g -^ Vg conversion acts in the opposite way. Depending on 
the 6m^ range, one channel dominates over the other one. For 5m? > 10^ eV^, only z/g — ^ i^s 
can occur which increases Y^ with respect to the case with no anti-neutrino conversion. For 
smaller values of 6m? there is an interplay of both conversions z/g — > i^s and u^ -^ Ug which 
can make Y^ < 0.4, hence enhancing the r-process. 

Properly averaging the neutrino absorption rates over the neutrino Fermi distribu- 
tion, we have calculated the electron abundance Y^. at the site where the heavy elements 
nucleosynthesis should take place as a function of {Sm,"^ , sir? 29) . In Fig. 7 we present our 
result. 

For a successful r-process, the region above Ye > 0.5 is ruled out. For Sm? > 10^ 
eV^ and sin^ 29 > 2 x 10^^ — 10~^ only the Ug — > z/^ channel undergoes the conversion. 
As we have discussed in Sec. 4, this region is also disfavoured by the neutrino re-heating 
consideration (Fig. 5). In correspondence of the corner delimited by 50 eV^ < Sm"^ < 10^ 
eV^ and sin^ 29 > 10^'^ the conversions take place in both channels z/g -^ Vg and z/g ^^ Vg 
and are adiabatic enough. In this case the significant suppression of the Ve flux and the 
recovery of the original z^e flux at r^ would lead to Y^ > 0.5 — 1. 

On the other hand we find that the supernova nucleosynthesis could be enhanced 
in the region enclosed by the dotted contour Y^ = 0.4, delimited by Smi'^ < 10^ eV^ 
and sin^ 26* > 10~^. Inside the contour of Ye = 0.33, Ye gets stabilised to 1/3 due to 
the feedback effect (in the absence of the feedback Ye would be lower than 1/3 - see the 
discussion in Sec. 3.2). We see that the mass range more promising for the neutrino hot 
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dark matter scenario, dm? < 10 eV^, is favourable for the r-process nucleosynthesis and it is 
neither in conflict with the re-heating process (see Fig. 5) nor with SN1987A observations 
(see Fig. 6). 



5 Discussion 

In this paper we have investigated the effect of resonant conversions of v^ or z/g into sterile 
neutrinos in the region above the hot proto-neutron star in type II supernova. For cos- 
mologically interesting mass values (1-100 eV) and mixing angle sin^ 2^ > 10~^ ^ 10~^, 
both z/e and z/g could be converted into Vs and Vg (respectively) in the region outside neu- 
trinosphere due to the non-monotonic behaviour of the effective matter potential. Such 
conversion could lead to the depletion of z/g and z/g fluxes, resulting in a suppression of the 
neutrino re-heating behind the stalled shock and of the expected z/g signal in terrestrial 
detectors. On the basis of these arguments we have derived constraints the neutrino mass 
and mixing parameters. We have found that for 5m^ sin^ 26 > 10~^eV^, the energy deposi- 
tion by z/g and z/g absorption reactions during the shock re-heating epoch (t < 1 s after the 
bounce) could be significantly decreased. This is not welcome for the delayed explosion 
scenario which relies on the revival of the shock by neutrino re-heating. 

The successful observation of the SN1987A z/g signal in the 1MB and Kamiokande 
detector can also be used in order to rule out a similar parameter range. Indeed, requiring 
that the total z/g flux during the thermal neutrino emission epoch (t ~ 1 — 10 s p.b.) 
should not be signiflcantly depleted by z/g -^ z/g conversion enables us to rule out the range 
5m^sm^2e> lO'^eV^ 

We further point out that depending on which conversion channel, z/g -^ z/g or z/g -^ z/g 
is dominant, r-process nucleosynthesis, which might take place in neutrino-heated super- 
nova ejects at 1 < t < 20 s p.b., could either be suppressed or enhanced. For the parameter 
range Sm^ sin^ 29 > lO^^eV^ where z/g -^ z/g conversion is dominant, Y^, at the nucleosyn- 
thesis site could become larger than 0.5 and hence r-process would be prevented, leading 
to the exclusion of this range of mass and mixing. On the other hand, r-process nucle- 
osynthesis could be enhanced due to the decrease of Y^ down to the minimum value 1/3. 
This is due to the fact that the z/g — * Ug conversion is dominant if the parameters are in 
the region 6171^ < 100 eV^ and sin^ 26 > 10"^. 

We have also discussed that, in contrast to the usual resonant conversion among 
active neutrinos (such as z/g ^^ z/^, z/^), the decrease or increase of Y^. due to z/g -^ Vg or 
^e -^ ^s conversion could be important in the estimation of the conversion probabilities. 
Indeed, the conversion could be suppressed when Y^ reaches a value close to 1/3. Such 
effect may take place both for z/g -^ z/g or z/g — * z/g resonant conversions. However, such 
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feedback effect should not be operative when the conversion occurs in the region where the 
electron or positron capture reaction is dominant over the neutrino absorption reaction, or 
if both conversions (z/g and Ve channels) occur in the same region. 

Finally, we wish to remark on the importance of the supernova constraints on active 
sterile neutrino conversions we have derived here. Most high-energy particle physics ex- 
periments are insensitive to the possible existence of sterile neutrinos since these do not 
couple to the electroweak currents. However, their admixture in the charged current weak 
interaction could show up in reactor neutrino disappearance searches. The laboratory lim- 
its on the mixing between electron neutrinos and sterile neutrinos are quite weak for this 
mass range. From reactor experiments the bound is weaker than sin^ 2^ < 0.01 PBI, far 
weaker than the supernova limits we derive. In contrast, for large mixing the reactor limit 
on 5m? < 10^^ eV^ is stronger than the corresponding supernova limits. 

A much stronger though also more uncertain argument to constrain active-sterile 
neutrino conversions comes from big-bang cosmological nucleosynthesis. Assuming that 



the number of effective neutrino species is bounded to be less than 4, one has p7[] : 

5m^ sin^2^ < 5 x 10"^eV^ , z/ = i/g, (1^ 



6m^ sin^20 < 3 x lO'^eV^ , z/ = i/^ , i/^ . (19) 

However the recent discrepant observational determinations of the primordial deuterium 
abundances [^, il] may force us to revise big-bang nucleosynthesis constraints on nonstan- 



dard neutrino physics. As a result the assumptions on which the previous limit is based are 
under debate [Q. It has been argued that Ny > Ai^ acceptable |3^, in which case there 



would be room to bring an extra sterile neutrino species into equilibrium with the known 
neutrinos in the early Universe and therefore no constraints on active to sterile neutrino 
oscillation parameters. Therefore, from this point of view, the restrictions on active-sterile 
neutrino oscillation parameters obtained from supernova theory and observations become 
quite relevant. 
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Fig. la: Typical matter density (solid line) and Fg (dotted line) profiles versus the radial 
distance from the center of the star, in Wilson's numerical supernova model at t = 0.15 
s after the core bounce. The diagonal cross indicates the position of the surface of the 
neutrinosphere. 
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Fig. lb: Same as in Fig. la but for t ~ 6 s after the core bounce. 
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Fig. 2a: The modulus of the matter potential Ve in eq. @ versus the radial distance r 
from the center of the star. This is obtained using the matter density and Yg profiles in Fig. 
la, at t < 1 s after the core bounce. The solid and dashed lines correspond to positive and 
negative potential, respectively, and the position where K = is denoted by r*. We also 
indicate, in the right ordinate, the Sm^ values for which a E = 10 MeV neutrino undergoes 
resonant conversion, for the corresponding value of \Ve\ on the left ordinate (small mixing 
angle is understood). 
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Fig. 3: Schematic figure illustrating the level crossings for v^ — ^s and v^ — ^s system. The 
energy levels of Ve (solid line), Ve (dashed line) and z/g, Vg (horizontal line) are given as 
functions of the stellar radius. The positions where Ve -^ ^s and Ve ^^ ^s resonances occur 
are indicated by ri, r2 and rs, respectively. 
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Fig. 4a: Main neutrino reaction rates (left ordinate) versus the radial distance from the 
stellar center at t < 1 s after the core bounce: z/g^ -^ pe~ (solid curve), v^P -^ ne^ (dot- 
dashed curve), e~p — > v^n (dashed curve) and e^n — > v^.p (dotted curve). The temperature 
profile (right ordinate) is also shown by the solid line (labelled T). We assume a neutrino 
(and anti-neutrino) luminosity of L^ = 10^^ ergs/s. 
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Fig. 4b: Same as in Fig. 4a but for t > 1 s after the core bounce. 
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Fig. 5: Contour plot of the ratio R of the neutrino energy deposition behind the shock wave 
in the presence of conversions into sterile neutrinos, versus the case without conversions, 
as defined in eq. (15). 
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Fig. 6: Contour plots of the survival probability P (figures at the curve) for the z/g -^ ^s 
conversion at t > 1 s p.b. The region to the right of the curves can be excluded by the 
observation of the SN1987A z/g signal. 
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Fig. 7: Contour plots for the electron concentration Y^ (figures at the curves) taking into 
account z/g -^ ^s and Ve -^ ^s conversions at t > 1 s p.b. The region to the right of 
the solid line labelled 0.5 is ruled out by the condition Y^ < 0.5 necessary for r-process 
nucleosynthesis to occur. For the parameter region inside the Yf. = 0.4 dotted contour 
r-process nucleosynthesis can be enhanced. 
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